ON THE RELATION BETWEEN DARBOUX 
TRANSFORMATIONS AND POLYNOMIAL MAPPINGS 



MAXIM DEREVYAGIN 



Abstract. Let dfi be a probability measure on [0, +00) such that its moments 
are finite. Then the Cauchy-Stieltjes transform S of dfj, is a Stieltjes function, 
which admits an expansion into a Stieltjes continued fraction. In the present 
paper, we consider a matrix interpretation of the unwrapping transformation 
S(A) I— > A5(A^), which is intimately related to the simplest case of polynomial 
mappings. More precisely, it is shown that this transformation is essentially 
a Darboux transformation of the underlying Jacobi matrix. Moreover, in this 
scheme, the Chihara construction of solutions to the Carlitz problem appears 
as a shifted Darboux transformation. 



1. Introduction 

It is well known that the monic Hermite polynomials iJ„ can be expressed in 
terms of the monic Laguerre polynomials Ln and Ln^'^'^ by means of the for- 
mulas '2T 

if2„(A) = 4-i/')(A2), iJ2„+i(A) = A4i/2)(A2), neZ+ = {0,1,2,...}. 



It turns out that this construction, which unwraps the measure, can be generalized 
to the case of monic polynomials Pj orthogonal with respect to a measure d/i on 
[0, +CX)) [6], [17]. In other words, we can complete the system {Pj{^^)}fLo to a 
system of monic orthogonal polynomials with the help of kernel polynomials Pj , 
i.e. the monic polynomials 



A 

which are orthogonal with respect to the measure tdfj,(t) on [0,+cx)). Namely, the 
sequence defined by the relation 

(1.1) $2,(A)=P,(A2), $2,+i(A) = AP,-(A2), jeZ+, 

is a sequence of orthogonal polynomials [5], [H]- Moreover, the measure of orthog- 
onality for these polynomials is the following measure 
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on R [9], [26]. This can be easily checked by recahing that the unwrapping trans- 
formation Pj I— >■ allows to reduce a Stieltjes moment problem to a Hamburger 
one |26j . |30| . Actually, the transformation in question can be expressed in terms 
of the Cauchy-Stieltjes transform 



Interestingly, iterations of the above simple construction can lead to orthogo- 
nal polynomials on several intervals. Moreover, fixed points of the unwrapping 
transformation give rise to orthogonal polynomials on Cantor sets and have non- 
trivial applications in physics [3]. Furthermore, the idea of unwrapping measures 
generated a study of polynomial mappings [I], [5], [IS], [21]. 

Concerning the completion problem (jl.ip . in [5], it is shown that it is not, in 
general, unique. More precisely, if the measure d/j, is supported on [a, +co) then 
there are infinitely many ways to complete the system {Pj{^^)}'j!Lo to a system of 
monic orthogonal polynomials. A construction of a family of completions to this 
problem was given in [9] by T. S. Chihara. The Chihara construction is the key 
ingredient to introduce and to study SDG maps |16j , |15j , which relate polynomials 
on the real line and the unit circle. 

The aim of this paper is twofold. The first goal is to give a matrix interpre- 
tation of the above constructions including the Chihara construction. The second 
purpose is to show that the generalized Jacobi matrices introduced and studied 
in [IT], [12], [in], [S] implicitly appear in the theory of polynomial mappings. An 
introduction to the theory of generalized Jacobi matrices from the classical point 
of view is also given. 

The paper is organized as follows. In Section 2, a preliminary information on 
the subclass N_oo of Nevanlinna functions is given. In particular, the subclass 
S-oo(c N_oo) of Stieltjes functions is also discussed. The goal of the next section 
is to introduce monic generalized Jacobi matrices associated with functions of the 
form F{X'^), where F G N_oo- By considering "symmetrized" generalized Jacobi 
matrices, Section 4 demonstrates the difference between the generalized Jacobi 
matrices associated with ^(A^) and with ^(A^), where F e N_oo and S G S_oo- 
In Section 5, it is shown that the unwrapping transformation S'(A) i— >■ \S(\'^), 
S € S_oo, can be realized as a Darboux transformation. Finally, this matrix 
interpretation is extended in Section 6 to include the Chihara construction [S]. 



In this section wc will give a brief description of the classical objects we are going 
to exploit in the present note. 

At first, recall that a Nevanlinna function is a function, holomorhic in the upper 
half-plane C+, which maps C+ into itself [1], [22]. We will denote the class of 
all Nevanlinna functions by N. Actually, in what follows we will only deal with 
Nevanlinna functions that admit the following asymptotic expansion 




of d/i as follows 




2. The Stieltjes class and the unwrapping 



(2.1) 




DARBOUX TRANSFORMATIONS AND POLYNOMIAL MAPPINGS 



3 



for all n G Z+ , where Sj , j G Z+ , are some real numbers (here A^oo means that A 
tends to oo nontangentially, i.e. inside the sector e < arg A < tt — e for some e > 0). 

As is known [T], [35], all Nevanlinna functions F e N_oo can be represented as 
follows 

where /i is a bounded measure supported on the real line JR. Moreover, we have 
that Sj — Jj^Pdj.i(t), j G Z+, are the moments of the measure d/Lt [I]. Throughout 
this paper, we assume that dfi is a nontrivial probability measure^ that is, sq = 1 
and det(sj+fe)",.^g > for all n G Z+. 

Once we have a measure with finite moments of all non-negative orders, it is 
possible to construct a sequence of orthogonal polynomials. It is well known that 
the sequence {Pj}^o of monic polynomials orthogonal with respect to jjL satisfies 
a three-term recurrence relation [T] , [3T] , [50] 

(2.2) AP,(A) - P,+i(A) + h,P,{\) + c,_iP,_i(A), J G Z+, 
with the initial conditions 

P_i(A)=0, Po(A) = l, 

where bj G M and Cj > 0, j G Z+. In fact, the relation (I2.2p can be rewritten as 
follows 

(2.3) Jp{\) = Ap(A), 

where p — (Pq, Pi, Pj, . . . )^ and J is a semi-infinite tridiagonal matrix of the form 

(b, 1 \ 

Co bi 1 

ci &2 ' ■ 

V •■/ 

The matrix J is said to be the monic Jacobi matrix associated with fi or, equivalently 
with F. It should be stressed here that F G N_oo uniquely determines J. However, 
J can correspond to infinitely many distinct functions of the class N_oo [1], [30) . 
In this case, by F we simply understand a family of N_oo-functions that generate 
the monic Jacobi matrix J. 

There is another way to get to Jacobi matrices and three term recurrence rela- 
tions. The starting point for this way is (|2.ip . More precisely, the asymptotic ex- 
pansion of P G N_oo can be expanded into the following continued fraction [2T], [30] 

1 1 Co Ci 



(2.4) J = 



(2.5) P(A) 



X-bo ^-^0 - A-6i 

A - 6i - — 



which is called a J-fraction. At this point, it is clear how to construct the Jacobi 
matrix J (|2.4p from the J-fraction (|2.5I) . Besides, the polynomials Pj appear in 
this scheme as denominator polynomials of the approximants to the J-fraction. 
Potentially, this method is more convenient for generalizations since we can start 
with a generic asymptotic expansion (which is not necessarily related to a posi- 
tive measure) and try to expand it into a good continued fraction (for instance, 
see [IT], [31], [33], [30]). Before going into details of such generalizations, let us 
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turn our attention to the Stieltjes functions which are basically the essence of the 
present study. 

We say that a function S € N_oo is a Stiehjes function if it has the following 
integral representation 

i.e. the support of the measure is contained in or equal to [0, oo). The class of all 
such functions is denoted by S_oo- 

Evidently, every S G S_oo has the following asymptotic expansion 

(2-7) '^^^-'i-^-----^+-{xi^)^ 
for all n G Z+ and the J-fraction expansion 

rV • 

A — Oo ^ A — oi — A — 02 ^ A — 03 — ... 

Remarkably, every Stieltjes function S G S„oo can be also expanded into the Stielt- 
jes fraction [21], [30] 

1 d^ dn d-i di di 

which is shortly called an S'-fraction. The fact that all functions of the class S_oo 
generate two different continued fractions, could be puzzling (this is not true for the 
entire class N_ooO- However, the things become transparent as soon as one recalls 
another equivalent definition of the class S-oo^ the function _F G N_oo belongs 
to S_co if and only if AF(A^) G N^oo [22]. So, there are two functions S{X) and 
AS'(A^) of the class N-oo related to every function S G S_oo and that is why S has 
two continued fraction expansions. Indeed, if S is of the form (|2.6p then one can 
easily check that 

(2.10) XS{X')^ j^-LJ!^d^,{n 

Therefore, one can write 

^^^^ ' A-/3o-A-/?i-A-/?2-A-/33----' 

Further, we see from (|2.10p that the corresponding measure is symmetric. This fact 
implies that /3j = 0, j G Z+ (see [1], [26], [30] for more details) and thus we have 

\c..\2N 1 70 71 72 73 

^^^^ ^^"A - T- T- T- T ■ 

Dividing the latter relation by A we get 

c./x2^ 1 70 71 72 73 

^^^^^ A^-T-A^-T-A^ ' 
which becomes (|2.9p after the substitution A^ — > A. Actually, the transformation 

(2.11) S{X) ^ A5(A2), 

which we will call the unwrapping by obvious reasons, is a way to reduce a Stieltjes 
moment problem to a Hamburger moment problem [301 Section 87] (see also [IHl 



2n 1 70 71 72 
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Theorem 2.13]). At the same tniie, this shows some importance of fmictions of the 
form S'(A^) in the theory of S'-fractions and J-fractions. 

3. MONIC GENERALIZED JACOBI MATRICES 

Now we are in a position to discuss another class of asymptotic relations and 
the corresponding matrices. Namely, we are going to associate a block tridiagonal 
matrix with a function ^ of the form 

(3.1) ;j(A)=F(A2), 

where F G N_oo- Hence, the function ^ admits the following asymptotic expansion 

oo 

(3.2) m=F{X')--Y.^, A^oo, 

j=o 

where ~ denotes the validity of (|2.ip for all n e Z+. Also, it will be useful to have 
all the powers in p.2p 

oo 

(3.3) 5(A) A^oo, 

where Sz] ~ and S2j+i = Sj for all j G Z+. As a matter of fact, the expansion on 
the right hand side of p.3p can be easily expanded into a P-fraction (see [21] for 
the general definition of P- fractions). 

Proposition 3.1. Let F be a function of the class N_oo- Then the function S^(A) = 
F{X^) admits the expansion into the following continued fraction (P-fraction) 

(3 4) — "1 

where Cj > and bj G M are the entries of the J- fraction (j2.5|) . 

Proof. The proof is straightforward from ()2.5p and p.l|) . □ 

Next, the basic theory of continued fractions ;21j reads that the denominators *Pj 
of the approximants to the fraction (j3.4l) satisfy the three-term recurrence relation 

(3.5) X^%{X) = V,+i{X) + b,V,{X) + c,_i*p,_i(A), J e Z+, 

with the initial conditions 

q3_i(A)=0, *Po(A) = l, 

where Cj > and bj E M. Obviously, deg^j = 2n and *Pj(A) = Pj(A^), where Pj 
are monic orthogonal polynomials generated by F. 

A natural desire coming from the classical setting is to represent (|3.5p as an 
eigenvalue problem for a matrix (i.e. to get a problem similar to (j2.3p ). The main 
obstacle to do that is the fact that p.Sp is not linear with respect to A. So, let us 
linearize the system (|3.5p . To this end, introduce a new family of polynomials: 

T2,(A)=*p,(A), l2,+i{X)^X^,{X), jeZ+. 

For the new sequence, (|3.5p takes the form 

AT2,(A)=T2,+i(A) 

Xl2j + l{X) =T2j+2(A) + bjl2jiX) + Cj_iT2j-2(A). 
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Now, the problem in question is linear and p.5|) can be rewritten as follows 






1 














... \ 






/^o(A)\ 


^'o 





1 













Ti(A) 




Ti(A) 











1 














^2(A) 


Co 





bi 





1 







2:3(A) 


= A 


^3(A) 

















1 




(A) 




T4(A) 








Cl 












(A) 




(A) 








■■• / 


V : ) 







V 

This brings us to the following concept. 

Definition 3.2 |14| . cf. |23j). A monic generalized Jacobi matrix associated 

with the function ^ is the tridiagonal block matrix 



/So 
£0 



(3.7) 



So 
£1 



Si 

232 



where the entries are as follows 

'0 1 







,2), = 



Actually, the monic generalized Jacobi matrix is constructed from the P-fraction 
p.4p . It turns out that one can associate a generalized Jacobi matrix of a similar 
type to any P-fraction (for more details see [TO l ITl ] [T2 ] ). 

Summing up we see that the P-fraction expansion of ^ leads to the following 
eigenvector problem 

at(A) = At(A), 

where i = (Xq, Ti, . . . )^. Keeping in mind the classical case, one would expect 
some orthogonality of {1j}°^o showing up. Mimicking the classical case, let us 
consider a linear functional 6 defined on the space V of all polynomials in A in the 
following way 



6(A^) 



where Sj are defined by 
with respect to ©, that is 



To begin with, observe that Tq is orthogonal to itself 



(3.8) 



©(To(A)To(A)) =©(1) =So = 0. 



This property demonstrates that © is not a regular functional. In other words, 
det(Si+j)f can vanish for some k g Z_|_. Besides p.8|) . the system {Tj}°^o 
possesses the following orthogonality property. 



Theorem 3.3. 

(3.9) 

(3.10) 



We have that 
©(X2,(A)X,(A)) 



D 



©(3:2,+i(A)Tfe(A)) = 



i-i 



j,k e Z+, 



where Dj = det(sijf^k)i k=o ^' -D-i = 1, Sfc are given by (j2.ip . and Sik is the 
Kronecker delta. 
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Proof. Recall that the monic orthogonal polynomials corresponding to F can be 
calculated by the formula 



P,(A) = 



1 



D 



So 



S2J-1 



1 Ai . 

We begin the proof of (13. 9|) and (|3.10|) by noticing that 



6(5:2, (A) A 



2l\ 



0, 



I e z+. 

The latter relation is true because T2j(A) = Pj{^^) (5(A^'') = S2r 
r = 0,1,. ..,21. Thus, we have just proved dSSJ for fc = 2? and (|3?T0)) for fc = 2/ + 1. 
The rest follows from the equality 



for 



6(T2,(A)A^'+i) 





So Sl 




1 








Sj_l Sj 


■ S2J-1 




S; S; + i . 





for/ = 0,l,...j. 



□ 



Evidently, one cannot say that the sequence {'Zj}^Q is orthogonal in the con- 
ventional sense. However, (I3.9p and (I3.10p show that the family {IjJj^Q is bi- 
orthogonal to itself with respect to &. M.G. Krein and H. Langer proposed to call 
such systems almost orthogonal in their paper about indefinite analogs of the Ham- 
burger and Stieltjes moment problems [23' and some earlier work. At the same 
time, the polynomials resemble skew orthogonal polynomials j2j and inherit 
many of their algebraic properties. 



Remark 3.1. From the proof of Theorem 

D 



(3.11) 



6(%(A)A'=) 
6(12,+! (A) A'=) 



D 



we can also get the following relations 
,2j + l, 



D 



-S2j + lk, fc = 0, 1, 

s. 



3-1 



0,1,. ..,2j. 



It should be mentioned that the first formula in p. lip is a special case of [25} 
Proposition 3.2]. 

Particularly, the results of this section are valid for functions of the form S'(A^). 
To analyze these functions we will need the following concept. 

Definition 3.4 ([13], [20]). Let us say that a function ^ is definitazable by the 
polynomial p{X) = A if it admits the asymptotic expansion 

S2n 



(3-12) A-+1 
with some Sj eR, j E Z+, and A^5(A) + Sq G N_o( 



A— ^00, 



Definitazibility is a very important concept in the theory of operators in Krein 
spaces [20], [24]. Generally, the concept means that a definitazable object, which 
acts in a Krein space, can be made definite (related to a Hilbert space) by means 
of some transformations. It should be also mentioned that the definizability by 
polynomials was also used in [T^] and in [3S]. 
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In view of (|2.10|) . the function 
(3.13) 3(A):=5(A2), 

where S G S_oo> is a definitizable function. It could seem from (|2.10p that 3 is 
the Cauchy transform of a finite signed measure. This can be easily disproved by 
considering the simplest Stieltjes function S'o(A) = — 1/A. Then 

3o(A)=5o(A2) = -^, 

which cannot be the Cauchy transform of a finite signed measure because it is the 
Cauchy transform of the first derivative (^j^j of the Dirac delta function. However, 
for example, if the support of the measure representing the original function S 
doesn't contain then 3 is the Cauchy transform of a signed measure. 

4. A BRIEF TOUR OF THE SELF-ADJOINT CASE 

This section demonstrates the difference between the generalized Jacobi matrices 
associated with ^ and with 3- To accomplish that, we have to start by considering 
the symmetric Jacobi matrix 

/ bo VcE \ 

Js 



corresponding to the function F . As a matter of fact, to realize the principal 
contrast between the matrices generated by ^ and by 3, it is enough to see it for 
bounded operators. So, in this section we can assume that Js is bounded in P . In 
other words, we restrict ourselves to the case of the measure c?/^ with a bounded 
support. 

Clearly, one can repeat the construction of Section 2 in the symmetric case and 
that gives the generalized Jacobi matrix 

\ 



/ 1 

6o 

















1 
6i 




V^r 












1 
52 




V '■ 









Evidently, Zs is not symmetric. However, the relations p.9p and (I3.10p suggest 
that we should consider the operator -Js in a different space. Namely, introduce the 
following inner product 

\x,y\ = {Gx,y), 

where (•, •) is the usual .^^-inner product and the Gram matrix G is defined as follows 

G = diag(Go,Gi,G2,...), = (^J , jeZ+. 

It should be stressed here that [•, •] is not positive definite: 

[eo, eo] = (Geo, eo) =0, eo = (1, 0, 0, 0, . . . )^. 
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Though it is not pleasant to work with indefinite inner products, the choice is 
dictated by the fact that we want the multiplication operator on V provided with 
the inner product defined by (3 to become the generalized Jacobi matrix Zs acting 
on the space P{G) of £^-vectors equipped with [x,y] = {Gx,y) (see [11], [23]). 
Notice the matrix G satisfies the following conditions 

G = G*, G^=I. 

Thus, the space £^(G) is a Krein space (for example, see [21]) ■ By definition, the 
norm of a Krein space is the same as the norm of the Hilbert space generating the 
Krein space [24] . In particular, the norm of the Krein space £'^{G) is the norm of 
the Hilbert space £'^. 

Proposition 4.1. The generalized Jacobi matrix associated with ^ defined 
by (13. ip is a bounded self-adjoint operator in i'^{G), that is 

[ZsX,y]^[x,Zsy], x,yee'^{G). 

Moreover, the spectrum a{Zs) of can be found by the formula 

(4.1) a{Zs) = {A e C : e a{J,)} cRUiR, 

where a{Js) denotes the spectrum of Js acting in the Hilbert space . 

Proof. The boundedness follows from the boundedness of J^., which we assume in 
this section. Next, one sees that the matrix 



bo 
1 











... \ 







bi 
1 


















62 
1 










•■ / 



is symmetric and this proves the self-adjointness. Finally, formula (14. ip is true 
because is a linearization of the problem {Js — \^I)x = y, x,y G i'^ (see also [lOl 
Theorem 5.3] for a more strict statement). □ 

Remark 4.2. For the bounded Jacobi matrix Js, it is well known that its m-function 
mj^(A) — {{Js — A/)^^eo,eo) coincides with F [I], [26], [30]. It appears that the 
same holds for bounded generalized Jacobi matrices [TU]. Namely, applying [TTl 
Proposition 2.8] to we arrive at 

"^0. (A) - [(aj - A/)-ieo, eo] - 5(A) = ^(A^), A G C \ a(3s). 

Actually, knowing that an operator is self-adjoint in a Krein space doesn't say 
much because the spectrum of a self-adjoint operator in a Krein space can be fairly 
arbitrary and the structure of the operator can be rather wild [24]. Therefore, it 
makes sense to study narrower classes of operators in Krein spaces. For example, 
one of such classes is the class of definitizable operators. In fact, the spectral 
calculus for definitizable operators is constructed j24| . 

Definition 4.3 ([24]). A bounded self-adjoint operator A in the Krein space £'^{G) 
is called definitizable if there exists a polynomial p such that 

[p{A)x,x]>0, xel^{G). 
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Now, we are in a position to formulate the following result. 



Theorem 4.4. The generalized Jacobi matrix associated with 3 given by (|3.13p 
is non-negative in i^{G), i.e. 



and, thus, 3s *s definitizable by the polynomial p{X) = 
that 

aids) = {XeC:\^e aiJs)} C 



A. Furthermore, we have 



Proof. It is easy to see that 



where = denotes the unitary equivalence in the Hilbert space £'^. Hence, the non- 
negativity of 3s follows from the non-negativity of Jg corresponding to the Stieltjes 
function (the statement about the non- negativity of Jg can be found in [26j). Then, 
it remains to apply (|4.1|) taking into account that Jg is non-negative. □ 



5. The unwrapping as a Darboux transformation 
The unwrapping transformation 
(5.1) S_oo 9 ^(A) ^ XS{X^) e N_oo 

plays an important role not only in mathematics |26| . [30) but also in physics [3]. 
It also gave rise to the study of polynomial mappings [7] , [5] , [H] , [13 . The goal of 
this section is to give a matrix interpretation of the transformation (jS.ip . Namely, 
we split (|5.ip into two steps: 

S_eo 3 Six) ^ 5(A2) = 3(A), 3(A) ^ XS{X^) e N_^. 

From the point of view of Jacobi matrices, the first step was considered in Section 2. 
Now we are going to interpret the second one. To be more precise, it will be shown 
that the second transformation is essentially a Darboux transformation (strictly 
speaking, a Christoffel transformation [29] . |31| ) of the underlying monic generalized 
Jacobi matrix obtained after the first step. 

We start by recalling the well-known facts about the Darboux transformations 
of monic Jacobi matrices [J]. Let S belong to S_oo, that is. 



SiX) 



d^i{t) 
t-X' 



Then the corresponding monic orthogonal polynomials obey 



(5.2) 

Indeed, we have that 



(^1)^P,(0) >0, jez+. 





So Sl . 




1 








Sj_i Sj . 






1 . 






D 



Sl 



S2J-1 
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and det{si+i+k)i'ii.^Q are positive since the moments Sfc ~ Sk+i, k e Z+, are the 
moments of the nontrivial positive measure tdfi{t) on [0, +cx)) [T], [33], [30] ■ Con- 
sequently, the corresponding Jacobi matrix admits an LJ7-decomposition T . This 
means that J can be represented as foUows 
(5.3) 



fbo 1 

Co bi 



Cl 



1 

62 



\ 



7 



LU 



/I 

h 



\ fuo 1 

Ml 







1 

"2 



V 



where the entries u„ and Z„ can be found by the formulas 

Pn+m 



(5.4) 



P„(0) 



> 0, In+i = — > 0, n e Z+. 



In fact, in this case, the factorization ()5.3|) is a version of the famous Cholesky 
decomposition, which is widely used in computational mathematics (for instance, 
see |18j). Moreover, it turns out that a similar statetment is also true for monic gen- 
eralized Jacobi matrices |141 Proposition 3.1] and the corresponding factorization is 
related to the Bunch-Kaufman decomposition [F . In the case of the monic gener- 
alized Jacobi matrices associated with functions of the form p.l3p , the statement 
is rather obvious up to the matrix structure. 

Proposition 5.1. Let ^ be a monic generalized Jacobi matrix associated with 3 
defined by (j3.13l) . Then T2j(0) = Pj{0) 7^ for all j G Z+ and -J admits the 
following factorization 

(5.5) a = 

where £ and it are block lower and upper triangular matrices having the forms 



(h 
£1 



(5.6) 




h 





h 



\ 



■J 



u = 



ill X>i 
il. 



in which the entries are of the form 



(5.7) 





u 





1 



I2 = 



1 
1 



Moreover, the following relations hold true 

„ _ ^2,+2(0) 



(5.8) 



12,(0) 



= Uj > 0, Ij+i ^ ^ ^ Ij+i > 0, j e 



where Uj and Ij are given by (|5.4p . 
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Proof. It is not so hard to see that the equahty 



/ ito £>o 



Si 

£2S)l +it2 



V 

reduces to (15.3 



/ V 



/'Bo 


So 




\ 




*Bi 


Si 






£1 


«2 




V 






■7 



□ 



Now, let us return to the monic Jacobi matrix J associated with the Stiehjes 
function S G S-a^- Recall [4], [29], [31] that the Darboux transformation of J is 
defined in the following manner 



(5.9) 



J = LU ^ J =UL. 



It is easy to check that the Darboux transformation J is a monic Jacobi matrix. 
Moreover, J is associated with the measure djl{t) ~ d^{t) [4] or, equivalently, with 
the function 

^(A) = A^(A) + so - XS{X) + 1= r 

Jo I — A 

Besides, from (j5.9p one can see how the corresponding orthogonal polynomials 
change under the Darboux transformation. Namely, in view of (15. 3p . the corre- 
sponding eigenvector problem (j2.3p can be rewritten as follows 

(5.10) UL{Up{X)) = X{Up{X)) 

and simple calculations show that 

P2(A) 



UpiX) 



^MX)\ 



P2{0) p 



iwy-^2(A) 



V 



Although the monic polynomials Pj{X) — p^^^^°Q^ Pj_i(A) satisfy the three term re- 
currence relation (|5.10p . they cannot be orthogonal since they all vanish at and, 
therefore, the sequence doesn't obey the proper initial condition. However, intro- 
ducing 

/h{x)\ 

Pi{X) _ 1 
^2(A) -A 

V ; / 

we get that 

JpiX) = XpiX) 

and Po{X) — 1. Consequently, Pj are the monic orthogonal polynomials corre- 
sponding to J — UL. Moreover, it can be easily seen from the orthogonality of 
Pj and the definition of the polynomials Pj , that Pj are orthogonal with respect 
to d'jj,{t) — td^{t). Actually, this is a special case of the Christoffcl formula ^27] 
Section 2.5]. 



m 



Up{X) 
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The Darboux transformations can be generalized to the case of generahzed Jacobi 
matrices associated with definitizable functions |14j . In particular, [14. Theorem 
4.2] can be complemented to the following statement. 

Theorem 5.2. Let^ be a monic generalized Jacobi matrix associated with 3 defined 
by 5W — S{X'^), 5e S_oo and let ^ = £,iX be its LU -factorization of the form (|5.5p . 
Then the matrix Z — is the monic classical Jacobi matrix associated with the 
function 3 given by 

3(A) = A3(A) +So = A3(A) = A5(A2). 

Besides, the orthogonal polynomials %j corresponding to Z can be calculated as 
follows 



(5.11) 



%(A) 

l2i+l(A) 



1 
A 



P. 



J el.. 



where Pj and Pj correspond to the monic Jacobi matrices J and J , respectively. 
Proof. First, simple computations show that 

\ 



/ 1 

Uo 




1 






... ^ 




( 1 ^ 

1 
















1 
ui 




1 





h 


1 
1 
















1 

U2 









h 


1 
1 




V '■■ 






■•• ) 


I ■■ 









is the tridiagonal matrix 



\ 



1 

h 



1 


Ui 



■I 



Taking into account (|5.8I) and Favard's theorem, we see that 3 is a monic Jacobi 
matrix associated with a positive measure. Next, let us see what happens with the 
corresponding polynomials. As in the classical case, it is also natural to introduce 
the following polynomials 



(5.12) 



t(A) 



5i(A) 

^2(A) 

V : } 



A 



ilt(A). 



Clearly, these polynomials obey the relations 

at(A) = At(A) 
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and To (A) = 1. Hence, 1j are monic orthogonal polynomials. It follows from (|5.12p 
that 

^2j(A) - |t2,+i(A), jeZ+, 

(5.13) ^ ^ 

T2,+i(A) = - (T2,+2(A) + w,T2,(A)) , J e Z+. 

Applying p. Ill) to (|5.13p . we get that are orthogonal with respect to the func- 
tional 6 defined in the following way 

6(/(A)) = 6(A/(A)), feV. 

Clearly, 6 generates the same moments as the function ^. This observation com- 
pletes the proof. □ 

We have just translated the unwrapping transformation 

^(A)h^3(A)^I3)(A) = A5(A2) 

into the language of matrices: 

(5.14) J = j(5)K^a = a(3)^a = 5(3). 

Basically, the essential part of (|5.14p is a Darboux transformation. So, we will call 
3 the extended Darboux transformation of J. 

It is worth mentioning here that 1-periodic dressing chains of the extended Dar- 
boux transformation of J, i.e. the equations of the form 

J-qI = Z, q>0, 

can lead to pure singular measures supported on Cantor sets These dressing 
chains are related to the Ferromagnetic Ising Model |2 . Some other dressing chains 
generated by Jacobi matrices can be found in [15], [28] . 
We conclude this section with the following. 

Remark 5.3. It is not so hard to see that an extended Darboux transformation can 
be also defined for the Jacobi matrix associated with F G N_oo provided that the 
corresponding monic orthogonal polynomials satisfy 

p,(o)^o, jez+. 

In this case, the matrix J = 3(A-F(A^)) is not necessarily generated by a posi- 
tive measure. Nevertheless, the matrix 3 can be constructed from the J-fraction 
expansion of AF(A^). 

6. The Chihara construction and shifted Darboux transformations 

Here we will give a matrix interpretation of the Chihara construction [9] (see 
also [15], [16], where it was used to study SDG maps) of solutions of the Carlitz 
problem related to polynomial mappings [B], |17] . For this purpose, consider the 
following function 

where a > 0. The gap (0,a) allows to do shifted Darboux transformations and, 
therefore, we can generalize Proposition EH] to the following statement. 
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Theorem 6.1. Let S have the representation (16. 1|) and let^ be a monic generalized 
Jacobi matrix associated with 3(A) ~ S'(A^), where S is given by (16. ip . Then for 
any a € (0, ^/a) we have that 

(~l)^T2,(a) = i-iyPj{a^) > 0, J e Z+, 

and ^ — al admits the following factorization 

(6.2) a - a/ = 

where £ and il are 6ZocA; lower and upper triangular matrices having the forms 



(6.3) 



/(So 
£i (£i 



V 



in which the entries are of the form 



(6.4) 





u 





1 



/Ho ®o 
ill 2)i 

il2 



-a 1 



1 

—a 1 



Moreover, the following relations hold true 



(6.5) 



P,(a2) 



> 0, [ 



■j+i 



> 0, i e 



Proof. Due to the structures of J, £, and il, the relation 

(a - a/) = £il 

is equivalent to the classical one 

which can be obtained from (|5.3p via considering the monic Jacobi matrix J^") = 
J — a^I corresponding to the Stietjes function 

dfj,{t) 



5(")(A) = 5(A + a2) = 



t-X 



(for more details see [4]). It remains to observe that 



D 



So 



Sj-i Sj 
1 a2 . 



S2J-1 



(-1)^- 



D 



si — a So 



a sq 



a'^Sj^i 



S2J-1 - a S2j-2 



and det(si+;+fc — a s;+fc); ^.^q are positive because Sk+i — a s^. are the moments of 
the nontrivial positive measure {t — a^)d^,{t) supported in [a, +cxi). □ 

The next step is to generalize Theorem 15.21 to the shifted case. 

Theorem 6.2. In the settings of Theorem \6.1\ we get that the matrix 

a^"' = il£ + a/ 

is the monic classical Jacobi matrix associated with the function ^'^'^^ given by 

^(")(A) = (A-a)5(A2). 
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Besides, the orthogonal polynomials xj"'* corresponding to 3'-"^ 
follows 

(6.6) 

5g|i(A) = (A + a)i>,(")(A2) 



be calculated 



A 



where the polynomials Pj correspond to the monic Jacobi matrix J associated with 
the Stieltjes function S . 

Proof. First, we can check by straightforward computations that the matrix 



v(") 



/ 


—a 1 














... \ 
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1 














... ^ 






Uo 


1 















—a 1 























—a 


1 















1 




















Ul 





1 









li 


—a 


1 










+al 













—a 


1 















1 




















U2 















(2 


—a 


1 






V : 








V : 






■•• ) 





coincides with the monic classical Jacobi matrix 



/-a 

Uo 



V 



1 

a 
h 



Ul 



1 

a 



■J 



Regarding the corresponding polynomials, we have the chain of transformations 
(a - aI)t{X) = (A - a)t(A) ^ (U£ + aI)m{X) = AUt(A) ^ a'"^ilt(A) = AUt(A). 
In this case, we have to normalize the vector lit in the following way 

t(")(a) 



t(A) 



in order the polynomials T^""* satisfy the proper initial condition 1f"'{X) — 1. The 
orthogonality with respect to the functional 

6(")(/(A)) = 6((A-a)/(A)), f eV, 
is also a simple consequence of p. lip . For example, let us check that 2^2^+1 (A) is 



ii-\x) 



V 



a 



-ilt(A) 



orthogonal to 1, A, 



A2^: 



G{xH^^Ux))=e(x 



^.+i(A^)-%§?^.(A^) 



= 6 (A^-p,+i(A2)) + u,6 (X^PjiX^)) - 0, fc - 0, . . . , 2i. 

Finally, the moments of 6^"^ coincide with the moments of (A — a)S{X^) given by 
the asymptotic expansion at infinity. □ 
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Remark 6.3. Interestingly, in Theoreni l6.1l and in Theorem l6.2l one can take ia, a € 
M, instead of a G (0, ^/a) and all conclusions will remain the same. In particular, 
we have that 



> 0, ij+i = ^ > 0, j e z+. 



Thus, the shifted extended Darboux transformation leads to the complex Jacobi 
matrix 





1 








ia 


1 






h 


~ia 


1 






Ul 


ia 



\ 



■J 



corresponding to the complex measure ^^^dfj,{t'^) on M (cf. [S]). 
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